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We compute the vacuum polarization on the lattice using non-perturbatively O(a) improved Wilson 
fermions. The result is compared with the operator product expansion (OPE). 



1. INTRODUCTION 

The vacuum polarization is given by 

IW?) = ijd 4 x e lf!X {0\TJ^x)J u (0)\0) 

= {QfiQu - q 2 g^)U(q 2 ), (1) 

with the vector current J^{x) — ijj(x)j^'ip(x). 
In the following we work at the Euclidean mo- 
menta Q 2 = —q 2 - Because the polarization 
function Q 2 ) has logarithmic divergences, 
it is customary to study the Adler function Q 
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D(Q 2 ) = -12tt 2 
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The Adler function provides a way of com- 
paring theoretical predictions from QCD with 
available time-like experimental data for the 
e + e~ total cross section via 



where 
R(s) = 
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is the ratio of the hadronic to the leptonic 
cross section at center of mass energy squared 
s. Lattice calculations will give valuable in- 
formation for the Adler function in the region 
shown in Figure |l|. 
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The standard operator product expansion 
(OPE) of the Adler function is of the form 

D(Q 2 )=D^ t {Q 2 )+D NP (Q 2 ). (5) 

The perturbative part is available to three 
loops in the form of a large-Q 2 expansion H H] 



and reads 
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To leading order in 1 /Q 2 the non-perturbative 
part is given by 
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The last sum in eq. (Q) is over q' dynamical 
quark flavors which will not concern us here 
since we are working in the quenched approx- 
imation. The Wilson coefficients are given by 
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Computation of the vacuum polarization al- 
lows us in principle to determine the strong 
coupling constant a s as well as the gluon, chi- 
ral and higher condensates. 

It has been claimed by several groups 
that the Adler function (as well as the polar- 
ization tensor) receives a further contribution 
of the form 



SD(Q 2 ) ~ A 2 /Q 2 



Such a term is not present in the OPE, because 
there exists no gauge invariant dimension-two 
operator. In this talk we shall compute the 
vacuum polarization and compare it with the 
predictions of the OPE. 

2. LATTICE EVALUATION 

To construct the polarization tensor on the 
lattice we demand the lattice Ward identity to 
be fulfilled. In lattice momentum space this 
leads to 



n M „ = 



(12) 



where the lattice momenta are defined as 
= (2/a)sm(q^a/2). (All momenta are Eu- 
clidean.) The polarization tensor splits into 
two parts: 



(13) 
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Figure 1. The analytic domain for the Adler 
function. Perturbative QCD is applicable in 
the hatched area of the complex q 2 -plane while 
lattice data can provide information in the 
shaded region. The picture is taken from ||/. 



The first part is given by 

H^ a \q) — e iqx+iq[ta/2-iq0a/2 
x 

x{0\J^(x)J^ + (0)\0), (14) 

while the second part corresponds to the tad- 
pole diagram 
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For the current in eq. ( |l4] ) we take the con- 
served vector current (CVC) 

4 1] ( x ) = + a V (1+7,0 U+(x)^(x) 

-$(x)(l-'y lt )U lt (x)1>(x + ap,)), (16) 



and the current in eq. (|15|) is given by 

4 2) (x) = + oA) (1 + 7,0 U+(x)i,(x)) 

+iP(x)(l- ll j)U l t(x)iP(x + a(i)). (17) 

We shall take the following ansatz for the po- 
larization tensor: 



(18) 



The calculations are done for non-pertur- 
batively O(a) improved Wilson fermions. Be- 
sides improving the fermionic action we have 
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Figure 2. The lattice data for (3 = 6.0, K = 
0.1345 on a 32 4 lattice. The line is the pertur- 
bative contribution. The scale is set using the 
r^-parameter pi 



Figure 3. The coefficient Pi of the power cor- 
rection 1 /Q 2 . The fit interval is chosen from 
approx. 1 GeV to 5 GeV for all lattices. 



to improve the current as well. The on-shell 
improved conserved vector current is 



+^-aid\ {'ip(x)a l j l \'ip(x)} 



(19) 



The value of the improvement coefficient c cvc 
is not known beyond tree level. We take the 
tree level value c cvc = 1. The derivative in 
eq. (pj|) must be defined such that the Ward 
identity, eq. is fulfilled. We take 

dxf(x) = ±{f(x + a\) - f{x - a\) + 
4a 

f(x + ajl + aX) - f(x + aft - aX)} . (20) 

In the momentum space this gives a factor 
—iq\COs(aq\/2)cos(aq^/2) which results in 
large 0(a 2 ) corrections, even in the free case. 
In the following we omit the 0(a 2 ) terms, thus 
leaving us with the factor — iq\. This keeps the 
Ward identity fulfilled. 



3. PRELIMINARY RESULTS 

We have performed simulations at the (3 and 
k values shown in Table [l| Because the Adler 
function involves a derivative with respect to 
Q 2 we prefer to work with the polarization 
function. The polarization function is thus 
written as 

-12tt 2 U(~Q 2 ) = c + a a s {fi 2 ) + c 2 a 2 s (fi 2 ) 



(21) 



We have included an additive constant (P4) to 
account for the logarithmically divergent con- 
tribution Q,[l(J- This depends on the renor- 
malization scheme. For the perturbative con- 
tribution we take the three-loop result |Q] 
renormalized in the MS scheme. The first co- 
efficient in eq. (^Tj) reads for a single quark 
with charge Q q = 1 
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Table 1 

The values of (3 and k used in the simulations. 
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and is related to the leading part in eq. (g) via 
a derivative with respect to ln(Q 2 ). For a s (/J. 2 ) 
we take the four-loop result jll],[l2| with 
Aj^g = 238(19)MeV The quark mass 

renormalizations are taken from Thus 
the perturbative part is completely known. 
We also allow for a 1/Q 2 contribution. Resid- 
ual 0(a 2 Q 2 ) corrections are accounted for by 
P 3 . In Figure |^ we show the lattice data for 
j3 = 6.0 and k = 0.1345. This is compared 
with the perturbative part of eq. (21) includ- 
ing a constant term P4. There we have also set 
fi = 1 1 a. However the result was found not to 
depend significantly on the exact choice of \i. 
For values of Q 2 > 2GeV 2 we find very good 
agreement between the lattice data and three- 
loop perturbation theory. Alternatively fitting 
eq. (f2l|), our data in lGeV < Q < 5GeV gives 
the result shown in Figure ||. It turns out 
that Pi is consistent with zero. (Note that 
r 2 , « 6/GeV 2 .) At present we are not able to 
quote a reliable number for P2. 
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The numerical calculations were performed 
on the Quadrics parallel computer of DESY 



